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THE  DYNAMICS  OP  SHELL 

ABSTRACT 

The  exterior  ballistic  motion  of  shell  is  considered  by  the 
methods  of  classical  mechanics.  These  methods  prove  to  be  powerful 
tools  for  a  qualitative  analysis  of  the  non-linear  motion  of  shell 
and,  In  certain  special  cases,  for  obtaining  quantitative  results 
with  a  minimum  of  computational  difficulty.  The  solution  thus  ob¬ 
tained  for  the  aircraft  gunfire  problem  has  proven  to  be  especially 
useful. 
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INTRODUCTION 


Interest  in  the  general  motion  of  shell  with  nonlinear  aero¬ 
dynamics  and  dynamics  has  been  academic  until  only  recently.  While 
the  first  attempt  at  a  study  of  the  general  motion  was  made  by  Fowler 
and  Lock  in  1922  (two  years  after  their  presentation  of  the  first 
definitive  work  on  the  subject  of  linear  motion'*'),  their  considerations 
showed  nothing  worth  intensive  study  since  they  only  considered  non¬ 
linearity  in  the  overturning  moment.  The  experimental  methods  were 
too  weak  to  disclose  any  good  data  on  the  Magnus  effect'  . 

4  5 

I.  L.  Synge  and  C.  H.  Kebby  have  made  similar  analyses  of  the 
shell  as  a  "top".  The  latter  reference  Is  extremely  detailed,  giving 
all  possible  cases  for  the  "Fowler  moment" . 

Not  only  was  the  experimental  evidence  lacking,  but  also  the 
practical  interest  in  large  yaw  motion.  For  most  weapons,  to  be  effect¬ 
ive  the  shell  must  be  kept  in  the  small  yaw  region.  This  is  a  problem 
for  linearized  theory.  Two  notable  exceptions  to  this  premise  are  high 
angle  artillery  fire  and  bomber  defense  gun  fire.  The  first  problem, 
popularized  by  the  requirements  of  jungle  and  hill  warfare,  would  most 
probably  be  best  solved  by  avoiding  large  sumraital  yaws,  but  this  does 
not  appear  possible.  The  second  problem  has  large  yaws  almost  by 
definition,  namely  the  launching  of  shell  into  crosswinds  whose  velocity 
is  comparable  to  the  muzzle  velocity.  The  second  problem  with  its 
urgency  has  supplied  the  necessary  funds  and  priorities  for  truly 
broaching  the  subject  of  large  yaw. 

In  1952,  L.  H.  Thomas^  considered  the  motion  of  shell  as  a  problem 
in  classical  perturbation  theory.  Unfortunately  the  heyday  of  classical 
techniques  was  at  the  turn  of  the  century,  and  few  if  any  present  day 
ballisticians  were  prepared  to  cope  with  the  mechanical  sophistication 
of  the  paper.  This  is  especially  unfortunate  since  the  classical  methods 
have  proven  to  be  so  powerful  In  handling  the  orbital  problems  of  the 
astronomer,  and  the  problems  of  orbits  with  their  multiple  periodicities 
are  not  unlike  the  problems  of  the  rigid  body  motion  of  shell. 
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» 


Fortunately  there  is  nothing  about  the  perturbation  methods  of 

classical  mechanics  which  requires  a  deep  understanding  of  complicated 

notions  of  physics.  This  was  stressed  by  Dr.  Thomas  in  conversation 

with  the  author.  That  is,  the  almost-constants  can  be  arrived  at  by 

either  considering  the  canonical  transformations  of  Hamiltonian  theory, 

the  conservation  of  energy  and  momentum^  or  just  ad  hoc  properties  of 

the  differential  equations.  In  fact  it  is  the  latter  approach  which 
7 

C.  H.  Murphy  applied  so  successfully  to  the  problem  at  hand. 

It  is  the  purpose  of  this  paper  to  reproduce  the  results  of  Thomas 
by  appealing  to  direct  considerations  of  the  equations  of  motion.  At 
the  same  time  it  is  desirable  to  correlate  the  direct  approach  with 
the  body  of  Hamiltonian  theory  so  that  the  numerous  methods  of  the 
astronomers  may  be  made  more  available  to  the  ballisticians .  Therefore, 
the  first  section  is  concerned  with  the  general  theory  of  canonical 
transformation  and  secular  variables .  The  remainder  of  the  paper 
develops  the  theory  of  the  motion  of  shell  independently  of  the  first 
section  but  with  identifying  references  to  the  first  section. 

A  departure  from  Thomas*  paper  is  the  use  of  the  function  defined 
by  equation  6.5  as  the  fundamental  motion  rather  than  a  sinusoidal 
approximation.  The  author  prefers  to  make  the  trigonometic  approxi¬ 
mations  as  late  as  possible  in  the  averaging  process.  Sections  12,  15, 
and  l4  are  extensions  of  Thomas’  paper. 

One  particular  limitation  of  this  paper  is  the  omission  of  the 
effects  of  gravity.  While  this  deletion  is  not  important  in  the  case 
of  flat  fire  and  rather  desirable  for  the  sake  of  clarity,  there  is  a 
large  area  of  interest  (Howitzer  fire)  in  which  the  effects  of  gravity 
are  of  prime  importance.  It  is  in  this  area  that  the  author  hopes  to 
extend  these  methods. 

1.  HAMILTONIAN  APPROXIMATION  THEORY- 

The  basis  for  classical  perturbation  theory,  and  for  that  matter 
any  approximation  theory,  is  the  study  of  the  variation  of  what  would 
be  constants  in  a  simpler  but  similar  system.  We  shall  consider  two 
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types  of  systems:  a  simple  rigid  body  with  a  conservative  moment  of 
force  and  additional  non-conservative  aerodynamic  forces  and  moments 
of  force  and  a  compound  shell  with  a  conservative  moment  of  force. 


* 


# 


For  the  first  type  we  assume  that  everything  possible  has  been  put 
into  Hamiltonian  form.  That  is 


P±  =  -  ^H/Sq.  +  f. 
CLj.  =  5h/Sp1  +  g± 


1.1 


where  the  q’ s  are  the  generalized  coordinates,  p’s  are  their  conjugate 

moments,  H  is  the  Hamiltonian,  and  the  f’s  and  g’s  represent  those 

* 

forces  and  moments  which  defy  a  Hamiltonian  formulation.  H  contains 
the  kinetic  energy  of  the  shell  and  those  aerodynamic  forces  which  are 
conservative  (usually  the  overturning  moment). 

We  further  assume  that  the  effect  of  the  f’s  and  g’s  is  to  produce 
a  small  secular  change  in  the  motion  which  is  typified  in  any  small 
region  of  time  by  H.  This  is  to  say,  the  conservative  motion  at  any 
instant  gives  the  frequencies. and  amplitudes  of  oscillations  in  the 
system  while  the  non- conservative  terms  produce  long  term  changes  in 
the  parameters  of  the  conservative  system.  To  isolate  these  parameters 
of  the  conservative  system,  we  shall  introduce  new  variables  such  that 
the  coordinates  are  almost  cyclic  and  hence  the  conjugate  momenta  are 
almost  constant.  These  new  momenta  will  be  the  desired  parameters. 


We  make  a  canonical  change  of  variables  as  follows: 


Let 


i  f  1 


1.2 


* 

For  a  general  discussion  of  the  Hamiltonian  formulation  see  References 
8  and  9*  These  are  only  two  of  many  possible  references. 
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We  shall  find,  that  for  our  purposes  the  case  with  all  but  one 
variable’s  being  cyclic  is  of  particular  interest.  That  is 
H  =  H(p^,  .  .  .  ,  p^,  q^).  To  complete  the  transformation  we  let 

pi  =  G(H’  P2’  '. *  *  >  Pn’  ql^ 
and  define  a  generating  function 

a. 


n 


,  Pq,  x)dx 


f  -  JZ  \  h +  f1  °<pi’  • 

i=2 

which  gives  the  transformation  equations 
Pi  = 

Q±  =  8F/8P. 


The  generating  function  gives  the  following  transformation  for 
the  new  coordinates: 

Qi  "  *i  +  f  9G/aPi  dx  i  f  1 

r1 

Q1  =  SG/SP1  dx  . 

Noting  that 

SH/aQ1  =  0 

Sh/BP^  =  8  ,  (the  Kronecker  delta), 

we  write  the  new  Hamiltonian  equations  as 
P1  =  (cSP^ap^fj  +  (dP^dqj)  Sj 

Q±  =  &i;L  +  ( BQ±/SPj )  f j  +  (dQj/dqj)  gj 
where  summation  is  to  be  taken  over  repeated  indicies. 


1.3 


1.4 


1.5 
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The  new  generalized  coordinates  are  the  time  in  the  conservative 
system  for  1=1  and  the  perturbations  produced  in  the  original  coordinates 
( <l)  ty  the  introduction  of  the  non-conservative  terms .  It  should  be 
pointed  out  that  no  approximations  have  been  made  up  to  this  point. 


Allowing  that  the  non-conservative  terms  may  be  considered  as  small 
and  hence  that  the  Q’s  for  i/l  may  be  taken  as  zero  and  -  t,  we 
obtain  the  secular  equations 


P.  =  h.  (P_,  .  .  .  ,  P  ,t)  1.6 

x  1  '  1’  nJ 

from  equations  1.5  by  substituting  the  solution  of  the  conservative 
system  in  the  right  hand  side.  A  further  simplification  can  be  effected 
by  replacing  the  right  hand  side  by  its  average  over  a  period  of  motion 
giving: 


t  +  T/2 

P1  -  (I/’S)  J  h.  (P^  .  .  .  ,Fn,  y)  dy  ,  1.7 

t  -  T/2 


the  P’s  being  assumed  constant  during  the  averaging  process.  T  is  the 
period  of  motion. 


Should  the  frequencies  of  motion  be  of  particular  interest,  a 
suitable  choice  of  variables  would  be  the  action-angle  variables.  If 
on  the  other  hand  the  magnitude  of  the  motion  is  of  first  importance, 
the  roots  of  the  energy  equation  would  be  more  appropriate.  The  latter 
will  be  the  case  for  most  of  the  further  discussion. 


The  systems  of  the  second  type,  mentioned  at  the  beginning  of  this 
section,  will  be  considered  in  the  last  section  since,  in  breaking  faith 
somewhat  with  the  introduction,  the  development  will  be  done  in  the 
Hamiltonian  manner. 


2.  THE  EQUATIONS  OF  MOTION 

As  was  mentioned  in  the  introduction,  the  following  development 
will  be  performed  in  a  direct  fashion,  i.e.  by  the  use  of  familiar 
equations  of  motion  and  simple  algebraic  operations  thereon.  Never¬ 
theless  the  development  will  parallel  the  Hamiltonian  method  outlined 
in  section  1,  and  at  appropriate  points  reference  will  be  made  to  the 
preceding  outline.  The  object  of  the  direct  development  is  twofold. 
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First  the  complicated  force  system  is  hard  to  identify  in  terms  of 
generalized  coordinates;  second  the  direct  approach  will  be  easier 
for  the  uninitiated  to  follow,  while,  at  the  same  time,  giving  a  view 
into  the  inner  workings  of  the  theory  if  continual  reference  is  made 
to  section  1  and  Textbooks  of  classical  mechanics. 

The  starting  place  will  be  the  equations  of  motion  in  almost 

the  same  form  as  presented  in  reference  10.  The  major  difference 

is  in  the  definitions  of  the  aerodynamic  forces  and  moments  for  which 

the  axial  component  of  the  velocity  is  used  in  the  reference  while 

the  total  speed  is  used  herein.  Also  we  shall  use  the  arc  length  of 

the  trajectory  for  the  independent  variable  rather  than  the  integral 

of  the  axial  component  of  the  velocity  .  The  effects  of  drag  and 
« 

spin  deceleration  'are  not  considered  although  they  may  be  added  with 
little  more  than  an  increase  in  bookkeeping.  The  effects  of  gravity 
are  not  considered,  which  limits  the  discussion  to  "flat  fire". 

Without  further  discussion  the  equations  of  motion  are  taken  to  be: 

X*  =  ip  i  -  l  JT  X 

L  2.1 

p»  =  iAvp/B  -  (V  JT  +  iJM)  k~2  X  -  JH  k2~2  p 

where 

X  Complex  orientation  of  the  velocity  vector  with  respect 
to  the  shell  axis 

p  Complex  transverse  angular  velocity  of  the  shell  in 
radians/caliber  of  travel 

*  d/ds 

s  arc  length  of  the  trajectory  in  calibers 

8  angle  of  yaw 

i  cos  8 

sin  8  |  X j 

A  axial  moment  of  inertia 

B  transverse  moment  of  inertia 

¥ 

Actually,  the  same  definitions  as  we  use  were  used  in  the  earlier  but 
probably  not  so  widespread  work  of  Kelley  and  McShane,  On  the  Motion 
of  a  Projectile  with  Small  or  Slowly  Changing  Yaw,  BRL  Report  446,  1944. 
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J 

p 

d 

m 

*L 

% 

Si 


axial  radius  of  gyration  in  calibers 

transverse  radius  of  gyration  in  calibers 

axial  spin  in  radians/caliber  of  travel 

pd^K/m 

density  of  air 

diameter  of  shell 

mass  of  shell 

lift  coefficient 

Magnus  moment  coefficient 

damping  in  pitch  coefficient 

static  moment  coefficient 


3.  THE  VARIABLES  H,  {,  and  J 


To  reach  the  starting  place  of  section  one,  we  must  change  the 
variables  of  section  two  into  canonical  variables.  Although  we  shall 
go  directly  to  a  set  of  variables  containing  H  (as  in  1.5),  it  is 
worth  while  indicating  the  intermediate  canonical  variables.  Using 
the  definition  of  Eulerian  angles  in  reference  8  with  5=9,  one  can 
obtain  a  set  of  canonical  variables.  The  of  section  one  would  be; 
Qg  =  fa  V±  =  P8J  P2  =  P2  =  p;  P5  =  Pj  =  Py  •  In  what  follows 


we  shall  use  the  dimensionless  §  and  v  in  place  of  P„  and.  P... 

d  3 

The  normalization  factor  is  Av.  In  this  section  a  bar  over  a  symbol 
represents  the  complex  conjugate. 


We  define  a  dimensionless  energy  (or  Hamiltonian) 

•  5-1 

which  is  the  sum  of  the  kinetic  energy  of  the  transverse  angular  motion 
and  the  potential  energy  of  the  overturning  (static)  moment 

v  ■  f  i  J-2 

where 

S  =  stability  factor  =  (Av)2 / (4b? ~)  • 

2 

The  normalization  factor  is  (Av)  /B, 
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We  further  define  a  dimensionless  angular  momentum  which  is  the 
component  of  the  shell’s  angular  momentum  about  the  velocity  vector 

§  -  217  (*£  +  M-0  +  5.3 

The  variables  H  and  JT  are  essentially  those  used  by  L.  H.  Thomas. 

Using  equations  2.1,  the  above  definitions,  and  some  straight¬ 
forward  algebra,  one  obtains  the  following  as  the  equations  of  motion 
in  the  new  variables: 

H1  =  "JHk2”2  (2H"2V)  -  JTki~2  (f  -  i)  +  Jg  JL  (!  -  *2) 

5.4 

I*  =  -(i  JL  +  JHk2"2)  (f  "  •«)  -  -  JL)  U  -  •g2) 

i1  *+  ^  ^(1  -  l2)  (2H-2V)  -  [§-Z)2\  JL  (1-/)  . 

The  radical  in  the  right-hand  member  of  the  last  equation  is  equal 

2 

to  G  of  the  first  section  multiplied  by  sin  5.  Of  course,  Z  is  not  a 
canonical  coordinate,  and  a  new  variable  ©,  which  is  more  closely  related 
to  the  coordinate  canonical  to  H,  will  be  introduced  in  section  6.  However 
Z  will  continually  occur  in  our  considerations,  and  a  word  about  the  use 
of  the  cosine  of  the  angle  of  yaw  in  the  definition  of  the  aerodynamic 
coefficients  is  in  order. 

It  is  common  practice  to  represent  the  aerodynamic  coefficients  as 
even  functions  of  the  sine  of  the  angle  of  yaw.  This  practice  is  satis¬ 
factory  for  angles  less  than  90° >  Beyond  90°  this  representation  gives 
the  value  of  the  coefficient  for  6  to  be  that  of  l8o°-  5  which  is  generally 
incorrect.  The  general  expression  for  an  even  function  (in  angle),  which 
is  not  completely  pathological,  is  a  Fourier  cosine  series  which,  if  the 
function  is  analytic,  can  be  represented  by  a  power  series  involving 
both  even  and  odd  powers  of  the  cosine.  Hot  only  does  the  cosine  appear 
as  a  logical  variable  in  the  dynamics  but  also  in  the  aerodynamics. 

Of  additional  interest  is  the  ease  with  which  a  cosine  series  can 
be  fitted  to  experimental  data  by  trignometric  interpolation^  and  further 
the  ease  in  converting  a  cosine  Fourier  series  to  a  cosine  power  series. 
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The  advantage  of  the  variables  H  and  JT  is  that  they  are  constant 
for  the  case  in  which  the  motion  is  that  of  a  top  with  a  generalized 
overturning  moment  and  no  other  forces.  However,  it  is  just  that  top 
motion  which  makes  the  equations  so  difficult  to  handle.  Essentially 
We  now  have  variables  with  most  of  the  high  frequency  top  motion 
stripped  out. 

4.  THE  FUNCTION  f(i) 

It  was  noted  in  the  last  section  that  the  radical  on  the  right 
of  the  last  of  equations  3 .4  was  the  function  G  of  section  one  mul¬ 
tiplied  by  the  square  of  the  sine  of  the  angle  of  yaw.  We  define  the 
function  f{Z)  as 

f(i)  =  (1  -  /)  (21  -  2V)  -  (J  -  Z)2  .  4.1 

This  function  dominates  much  of  the  motion  of  a  shell  in  the  same 
manner  as  the  similar  function  dominates  the  theory  of  Abelian  integrals. 
In  fact  if  V  is  approximated  by  a  truncated  cosine  series,  f  is  a  poly¬ 
nomial  in  Z,  and  the  analysis  depends  strongly  on  the  theory  of  elliptic 
and  hyperelliptic  integrals.  Without  the  use  of  such  sophisticated 
mathematics,  much  can  still  be  learned  about  the  overall  motion  from  the 
roots  of  this  polynomial. 

Since  the  cosine  of  the  angle  of  yaw  must  be  a  real  quantity,  it 
follows  from  (3*5)  that  the  function  f  must  be  positive  for  any  value 
of  yaw  the  shell  can  assume  consistently  with  the  initial  conditions. 

In  particular,  the  yaw  will  oscillate  between  two  neighboring  roots 
Zq  and  i  such  that  the  function  is  positive  between  these  two  roots. 
There  is  of  course  the  requirement  that  the  cosine  of  a  real  angle 
must  be  of  not  more  than  unit  magnitude.  In  this  connection  it  is  grat¬ 
ifying  that  the  function  is  non-positive  for  Z  of  unit  magitude. 

The  special  case  of  V  linear  in  Z  corresponds  to  the  common  grav¬ 
itational  top.  Tliis  case  is  discussed  at  length  in  reference  8  and,  for 
that  matter,  most  any  text  of  classical  mechanics. 
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The  roots,  £q  and  £^,  of  f  which  bound  the  motion  are  often  well 
adapted  to  use  as  secular  (slowly  varying)  variables.  Indeed  the  entire 
concept  of  stability  hangs  on  these  roots.  In  many  cases  it  might  be 
desirable  to  avoid  the  algebraic  difficulties  of  finding  roots  of  high 
order  polynomial  by  carrying  both  the  roots  and  H  and  JT  as  dependent 
variables. 

5.  THE  ROOTS  AS  VARIABLES 

In  this  section  we  shall  do  two  things,  find  the  differential 
equations  for  the  roots  and  evaluate  H  and  _Jf  in  terms  of  the  roots. 

The  first  is  necessary  for  using  the  roots  as  dependent  variables  and 
for  determining  the  behavior  of  the  variable  9  which  will  be  introduced 
in  the  next  section.  The  second  is  necessary  if  the  roots  are  to  be 
used  exclusively  as  dependent  variables. 

We  proceed  as  follows; 

Let  r  be  one  of  the  roots  of  interest.  Then 
f(r)  =  (l-r2)-(2H-V(r))  -  (J-rf  =  0 

Differentiating  this  expression  we  obtain  , 

(df(r)/dr)r*  =  2($-r)J*  -  2(l-r2)Hl 

=  2(J-r)  [-(i  JL  +  jRk2  )^)_(j^'-2_jl)(1_/ 

-2(l-r2)  [-JHV2  (2H-2V(i))  -  +  JL(l-/)/4S(i)]. 

Using  f(r)  =  0  and  observing  care  in  keeping  £  and  r  separate,  one 
obtains,  by  eliminating  terms  in  J^2, 

(df(r)/dr)r*  =  2JL  (r-i)+[2H-2V(r)]  (r-i)(l+ri) 

+  [ l/4S(r)  -  1/kS  (4)]  (1-/)  +  \  (1-/) 

+  2k1'2  JT  (r-je)|l+ri  -  J(r+i)j  5.1 

+  2k2"2  JH|(r-i)(r-i)  +  2(l-r2)  [v(r)  -  V(i) 

Each  J  is  a  function  of  £. 
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It  is  certainly  reasonable  to  assume  that  V  and  S  are  different¬ 
iable.  In  which  case,  using  the  mean  value  theorem,  we  can  write: 

(3f(r)/Sr)r»  =  (l-i2)  JL  (5f  (r)/c>r)  +  (r-i)C  (r,i)  , 

Where 

lim.  C(r,i)  is  finite. 

Z  — *-r 

Anticipating  the  next  section,  we  define 


f(i)  =  (Jri)(i-i0)  F(i) 


5.2 


where  F  is  bonded,  positive,  and  non- vanishing  in  the  closed  interval 
Uq,  ^1)‘  Then 


(df(JQ)/dr)  =  (iri0)  F  (iQ) 

(Sf^/dr)  =  (iQ-ii)  F  (iL) 
and  £  _£ 

A,  =  (1-i2)  JL  +  ^  c  <*0>h/*o 

h  ■  h--*2)  ^  C<*L>«^X 


where  F~  and  F,  are  the  values  of  F  evaluated  at  &  and  A . 
-01  o  1 


5-5 


A  more  detailed  evaluation  of  the  derivatives  of  the  roots  would 
depend  on  a  knowledge  of  the  functional  form  of  V.  For  our  purposes 
equations  5*5  are  sufficient. 

The  next  step,  as  indicated  at  the  beginning  of  this  section,  is  to. 
find  the  functional  relationship  between  the  roots  and  H  and  J$".  To  do 
this  most  simply  and  to  give  the  most  understandable  form,  we  introduce 
two  new  intermediate  variables.  These  new  variables  are  quite  closely 
related  to  rather  familiar  concepts  in  exterior  ballistics.  In  particular 
let  5^  and  5^  be  the  angles  of  yaw  associated  with  Z^  and  Z^.  We  then 
define  two  new  variables 


p  =  (1/2)  (80  +  81) 
n  =  (1/2)  (8Q  -  51) 


5.4 


15 


The  letters  p  and  n  are  chosen  for  the  connotation  of  precession  and 
nutation.  Here  we  are  at  variance  with  the  majority  of  the  literature 
in  the  field  of  exterior  ballistics,  but  in  keeping  with  the  traditions 
of  classical  mechanics.  That  is  in  this  paper  we  call  nutation  the 
motion  involving  changes  in  the  angle  of  yaw  and  precession  the 
average  motion  of  the  shell  axis  about  the  velocity  vector.  In  this 
sense  we  shall  see  that  nutation  is  always  more  nearly  associated  with 
the  ballistic  concept  of  fast  motion  while  the  precession  can  be  either 
fast  or  slow. 


From  the  two  equations 

f(iQ)  =  0  and  f(j  )  =  0 


eliminate  H  to  obtain 


§ 


i  +  Vl  .  J. 

a  .  />  ^  /]  .  n  V  - 


^ 


*0  +  *l 


1-4 


VJ1  w 

2 


where  VQ  and  are  the  values  of  V  for  iQ  and  i^. 

Changing  variables  to  p  and  n  gives 
cos  n  (l  +  cr)  cos  p  (l  —  q) 

J-- 


2  cos  p 


2  cos  n 


where 


a  =  V 1  -  4  cos  p  cos  n  (V1-V0)/(2.1  -  &0)'  , 


5.5 


5.6 


which  reduces  to  the  usual  q  in  ballistics  in  those  cases  where  V  is 
linear  in  the  cosine  (K^  is  constant).  Now  if  0  is  the  azimuth  angle 
as  used  in  Reference  8, 

0.  =  Av  ifzll 

Bsin2  8 


If  we  consider  the  case  of  pure  precession  (n  =  0) 

Av  (l  +  q) 
r  2B  cos  p 

For  the  case  of  statically  unstable  shell  (0<q<l),  there  are 
two  precession  rates,  both  in  the  same  direction,  one  fast  and  one 
slow.  For  the  case  of  statically  stable  shell  (q>l)  there  are  also 


5.8 
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two  rates  tut  now  of  different  sign  as  well  as  magnitude.  For  the 
case  of  neutrally  stable  shell  (a  =  l)  there  is  only  one  rate  and 
that  is  the  vacuum  rate  as  one  would  expect. 

One  final  observation  about  the  rates  is  of  interest.  For  p  =  9°°, 
there  is  only  one  possible  precession  rate.  This  motion  is  possible 
since  the  limit  of  a  as  p  approaches  90°  is  one.  The  slow  motion  is 
the  possible  one  and  in  this  case 

I  =  cos2n  (V1  -  VQ)/(i1  -  iQ) 

In  a  similar  manner  as  for  equation  5*5  we  obtain: 

2  2 

if  p  1  +  o  o  1  ^  ®  1  1 

H  =  |  |  tan  p  (— g— )  +  tan  n  {— g— )  j  +  g  (v1+Vo)  5'9 

6.  THE  UNIFORMISING  VARIABLE  0 

In  Section  5  it  is  noted  that  l  is  not  a  canonical  coordinate. 

While  time  (actually  arc  length)  is  the  coordinate  conjugate  to  the 
Hamiltonian,  it  is  more  convenient  to  introduce  another  variable  0  by 
the  relation 

£  =  (l/2)(i0+i1)  +  (l^K^-^)  cos  ©  .  6.1 


Thomas  calls  this  variable  a  uniformising  variable  since,  if  SL 
varies  sinusoidally,  ©  varies  linearly.  In  any  event,  the  transfor¬ 
mation  absorbs  much  of  the  oscillation. 


Differentiating  6.1,  substituting  in  3.4,  and  rearranging  terms  gives 

Av  Vm  -  (h2»s-9)+i;  (k£°2_S) 


^0  ~  ^1  * 

-  2  Sin  QQ  =±I 


1 


Using  5*2  and  5-3  we  obtain 

0  B  2  Uq-^7 


I- C(iQ,i)  £)' 


F 

0 


1 


6.2 


In  general  the  presence  of  the  term  in  the  denominator  of 

the  last  term  causes  trouble  when  the  roots  are  almost  equal.  A  dis¬ 
cussion  of  this  case  will  be  given  in  Section  8.  Fortunately  for  roots 
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which  are  separated  sufficiently  that  the  manner  of  finding  the  average 
of  l  is  important,  the  entire  term  can  be  disregarded.  This  latter 
simplification  is  applicable  only  for  the  case  of  heavi'er-than-air  pro¬ 
jectiles  (see  below).  Since  this  is  the  case  with  which  we  are  con¬ 
cerned,  we  shall  use 

9*  =  ~  Vf (if  6.3 

as  the  equation  of  the  fundamental  motion.  This  is  similar  to  the  use 
of  the  approximation  =  t  in  Section  1. 


The  frequency  of  nutation  can  be  defined  as  the 
time  (arc  length)  in  which  0  changes  by  2it.  We  call 
length)  the  period,  P,  where 


-r 


e  +2  it 
d  e 


reciprocal  of  the 
this  time  (arc 


6.4 


To  establish  under  what  conditions  equation  6.3  is  valid,  we  shall 
evaluate  the  C  and  F  terms  for  =  £^  =  £  =  cos  p.  We  shall  then  con¬ 
sider  under  what  conditions  the  term  Av  V  F(cos  p)/B  dominates  the  term 
C(cos  p,cos  p)/F(cos  p)  -  ij  . 

Wow,  using  5-5;  and  5-9;  we  obtain 

C(cos  p,  cos  p)  =  2J^  £tan^p  (l+ct)^/4  -  sin^p  (l-cr2)/4j 


-2  2 

+  2Jt  kx  sin  p  cr 


6.5 


_2 

-  2J^  kg  tan  p  sin  p  (l+cr)  / 2 

Equation  6.5  may  be  used  in  detail  to  obtain  an  estimate  of  the  size  of 
C.  However,  in  the  cases  which  are  usually  of  interest,  we  can  merely 
say  that  C  is  of  the  order  of  magnitude  of  (pd  /ra)sin'p,  and  hence 

is  of  the  order  of  magnitude  of  (pd^/m)(sin  p)/(sin  n).  The 

comparison  is  then  between  Av-/f  /B  and  (pd^/m) (sin  p)/(ein  n)F.  That 
is,  the  approximation  6.3  is  invalid  for  n  such  that 


sin  n<(pd^/m)  sin  p  /  [(A  /B)  ^  ]  . 


6.6 
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=  1  to  get  as 


For  spin  stabilized  shell  we  may  use  F  ~  1  and  K^k^ 
an  order  of  magnitude  approximation 

sin  n  < (l/bS) (Av/b)  sin  p 
for  the. approximation  6. 3  to  be  invalid. 


6.7 


A  more  specific  analysis  along  the  same  lines  is  necessary  for  any 
particular  fin-stabilized  shell, 

7.  THE  FUNCTION  F 

The  function  F,  together  with  Av/B,  determines  the  frequency  of 'the 
nutation  of  the  shell  (6.1i).  If  V  is  at  most  quadratic  in  Z,  the  nutation, 
can'be  expressed  in  terms  of  elliptic  integrals.  If.V  is  of  higher  order.' 
the  solution  requires  hyper- elliptic  functions  (or  series- expansions) l 
It  is  Instructive  to  reduce  F  to  a  form  similar  to  the  expressions  for 
the  nutational  frequency  found. in  the  usual  ballistic  theory,  (in  view, 
of  our  definition  of  nutation,  the  nutational  frequency  we  shall  obtain 
will  be  equal  to  the  difference  of  the  two  frequencies  of  the  usual 
ballistic  theory). 

We  first  note  that  if 

f(i0)  =  f(V  =  o 


Then,  by  using  the  mean  value  theorem. 


f (i)  =  -  (1/2 )(i0-i)U~O  ^4^ 


7.1 


<1Z 


where  Z  J.ies  between  ZQ  and  Z^  and  depends  on  Z. 
Letting 


we  have 


V  =  V(z)  etc. 
F  =  1  +  (2H-2V) 


dV  „2v 

^  d Z  +  ^ 


d2V 
d  / 


7.2 


Several  cases  are  simple  enough  to  be  of  interest.  First  consider  the 

case  of  small  nutation.  Letting  n  go  to  zero  and  using  7.2,  5.9,  and  5-6, 

we  get  the  following  expression  for  small  nutation. 

2 

q  q  1  "i  CT  q  .  2 

F  =  o'  +  tan  p  ( — * — )  +  sin  p  — ^  7.5 

d  Z 
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The  other  special  case  we  shall  consider  is  that  for  V  quadratic 

2 

in  Z.  In  particular  let  V  =  a.Z  +  bj8  .  Then 

2  2 
2  2  1  +  2  1  + 

F  =  cr  4-  tan  p  ( — ■= — )  +  tan  n  ( — ~ — ) 


1  -  a 


(1  - 


2 
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7.4 


zl  +  *o 


)  +  2(1  -  Z  )  b 


8.  THE  BEHAVIOR  FOR  SMALL  NUTATION 


In  developing  the  equations  for  the  uniformising  variable,  mention 
was  made  that  the  right  hand  side  of  equation  6.2  could  cause  trouble. 
Such  is  often  the  case  when  the  nutation  is  quite  small  (while  a  specific 
check  would  be  required  for  each  case,  it  appears  that  generally  the 
approximation  presented  by  6.3  is  valid  for  any  nutation  which  can  be 
considered  sensibly  different  from  zero.  This  appears  to  be  the  case 
as  long  as  the  density  of  the  projectile  is  large  compared  to  the  density 
of  the  resisting  medium.)  While  the  case  of  extremely  small  nutation 
presents  no  problem  as  to  the  choice  of  a  proper  averaging  technique,  for 
completeness  it  is  desirable  to  discuss  this  case.  Indeed,  in  one  sense 
at  least,  it  is  mandatory  that  an  explanation  be  given  as  to  why  a  per¬ 
turbation  technique  appears  to  break  down  under  those  conditions  for 
which  the  perturbed  variables  should  be  changing  the  least.  It  is  de¬ 
sirable  to  show  that  the  difficulties  are  only  in  the  coordinate  system. 

Actually  the  variables  H  and  do  change  slowly  at  all  times.  The 
ficulty  is  in  the  roots  of  f.  When  the  two  roots  are  close  together, 
very  small  changes  in  the  parameters  of  the  polynomial  cause  proportion¬ 
ally  much  larger  changes  in  the  roots.  What  happens,  as  we  see  below,  is 
that  in  certain  cases  the  uniformising  variable  has  limited  variation, 
and  the  roots  oscillate  markedly. 


For  the  case  of  small  nutation  equations  5*5  and  6.2  can  be  approxi¬ 
mated  by 

-  -  ^ 

8.1 

=  cos  G  c/f 


where  C,  F,  and  q  are  evaluated  for  =  JJ  . 
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These  equations  can  he  normalized  by  introducing  the  variables 


x  = 


dz  = 


(i0  -  \)  CAv  7F 


FB 

Av  y^F  ds 
B 


The  normalized  equations  are 


8.2 


dO  sin  Q 

dz  x 


dx 

dz 


cos  0 


8.3 


A  graphical  representation  of  these  equations  appears  in  Figure  1. 

It  can  be  seen  that  two  types  of  solutions  are  possible:  closed  paths 

1  It 

around  the  points  x  =  1,  0  =  (2k  +  — )  or  around  the  points  x  =  -1, 

©  =  (2k  +  j  and  open  paths  in  which  ©  increases  indefinitely.  Actually 
the  variables  x  and  ©  are  not  the  observable  quantities.  What  is 
observable  is  the  change  in  the  .cosine  of  the  yaw,  this  change  being 
proportional  to  x  cos  ©.  If  we  denote  this  quantity  by  y  then  8,1  gives 


y  =  0 


which  by  the  definition  of  z  means  that  the  cosine  of  the  yaw  has  the 
usual  frequency,  viz.  Av 

9.  THE  SECULAR  EQUATIONS 

While  equations  6.3  and  the  first  two  of  equations  3. t  could  he  used 

their  present  form  for  machine  computation,  the  presence  of  high  frequency 

oscillations  is  a  source  of  possible  divergence  of  the  truncation  error 

in  the  computational  techniques  unless  extremely  small  integration  inter- 
12 

vals  are  employed  .  Further  the  presence  of  three  variables  upon  which 
the  equations  depend  explicitly,  rules  out  the  possibility  of  a  geometric 
interpretation  with  any  intuitive  value.  A  simplification  of  the  equations 
is  therefore  indicated  from  both  the  analytic  and  computational  views. 


21 


To  this  effect  we  invoke  the  sane  approximation  which  was  used  in 
obtaining  equation  6.3  from  6.2,  namely,  the  transient  times  are  large 
compared  to  the  nutational  period.  The  secular  equations  (see  1.6)  are 
obtained  by  averaging  the  dependence  of  the  first  two  of  equations  3.4 
on  0  over  a  period  of  nutation.  It  is  assumed  that  over  this  period  H 
and  jf  are  constant.  The  average  will  depend  upon  the  functional  form 
and  roots  (i^  and  of  f. 

Somewhat  more  specifically  the  modus  operandi  is  as  follows: 


Equations  3.4  may  be  written  as 

H’  ‘  all  H  +  a12  J+  a 

jT’  =  a00  §  + 


13 


9.1 


22  23 

where  the  a’s  are  functions  of  Z  and  hence  of  0  and  the  roots  of  f.  The 

a’s  are  replaced  hy  their  averages 

s+P 


J  a^0’h’  °)d£ 


9.2 


i 

p 


2n 


a {ZQ>  Zv  0)  d© 


Av 

B 


-/F(i0,  Z~9.1 


2m 

p.  f  -  39  -  _ 

The  roots  (i^  and  jg  )  are  considered  constant  in  these  integrations. 
Usually  the  a’s  can  be  written  as  polynomials  in  Z  and  hence  as 


hj  (1)  =  suk ( Vl  oosli  9 

with  summation  over  k. 


9-3 
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Therefore 


where 


ijk  bk 
cos^Q  cLQ 

~7? 

2jt 


9-4 


Evaluation  of  the  latter  quantity  involves  complete  elliptic  integrals 
of  the  first  and  second  kind  when  f  is  cubic  (see  next  section),  complete 
elliptic  integrals  of  the  first,  second,  and  third  kind  when  f  is  quartic 
and  complete  hyperelliptic  integrals  when  f  is  of  higher  order. 

The  final  secular  equations  are 


H’  ~  allkbkH  +  &  a12kbk  +  ai3k\ 

9-5 

=  a22kbl/  +  a23kbk 

where  the  a’s  and  b's  are  functions  of  ZQ  and  Z^. 

We  may  now  think  of  the  dynamic  history  of  the  shell  as  being  a  path 
in  the  H,  jH  plane.  A  shell  may  be  defined  as  dynamically  stable  in  an 
H,  region  if  all  paths  from  that  region  lead  to  the  H,  JT  point  corres¬ 
ponding  to  Zq  =  J&i  =  1 .  The  region  of  physically  consistent  H  and  § 
(that  is  H  and  §  such  that  f  is  positive  for  some  values  of  Z  in  the 
interval  -1<Z  <•  1 )  is  bounded  by  the  curve  of  values  of  H  and  j5"  for  which 
the  roots  are  equal  (n=0).  Any  limit  points  on  this  curve  represent 
stable  pure  precessional  motion.  Limit  points  inside  the  region  rep¬ 
resent  stable  "elicyclic"  motion.  In  those  cases  where  the  curves  for 
H*  =0  and  for  JT1  =  0  can  be  obtained,  the  intersections  of  these  curves 
(singular  points)  may  be  analysed  and  an  overall  description  of  the 
motion  obtained  (in  the  same  manner  as  reference  7).  The  latter  situ¬ 
ation  is  not  too  likely  since,  even  if  the  elliptic  integrals  can  be 
approximated,  a  system  of  three  nonlinear  algebraic  equations  must  be 
considered  (two  for  the  vanishing  of  the  derivatives  and  one  for  the 


roots  of  f).  A  more  profitable  line  of  attack  would  probably  be 
to  lay  down  a  sufficiently  dense  set  of  numerical  solutions  for 
reliable  inferences  to  be  drawn  about  the  general  motion. 


However,  there  are  many  special  cases  in  which  such  things  as 
restricted  initial  conditions,  etc.  allow  analytic  approximations  to  be 
made  and  a  same what  closed  solution  to  be  obtained. 

10.  THE  CASE  OF  A  CUBIC  f 

For  the  case  in  which  the  stability  factor  (S)  is  a  constant, 

(V  is  linear  in  &  as  in  the  usual  gravitational  top),  the  function 
f  is  a  cubic  polynomial  in  &,  and  the  secular  equations  can  be  written 
in  terms  of  complete  elliptic  integrals  of  the  first  and  second  kind. 

In  fact,  only  one  transcendental  function  need  be  used,  namely  E/K. 

To  demonstrate  this  fact  and  to  derive  the  detailed  equations  for 
the  cubic  case,  it  is  convenient  to  use  all  three  of  equations  3-4. 
Averaging  with  respect  to  2  is  equivalent  to  averaging  with  respect 
to  the  uniforraising  variable  9. 


where  again  the  roots  of  f  are  assumed  to  remain  constant  during  the 
integration. 

Wow 

f  =  (1  -  i2)(2H  -  1/23)  -  (J-£)2 

=  (1/2S)  (2  -  £0)(£  -  2±)(2  -  i2)  10>s 

t2  =  2S(2H  +  1)  -  tQ  -  2± 

Noting  that  f(l)  and  f(-l)  are  negative  and  that  since  f  is  an  odd 
order  polynomial  It  must  approach  plus  infinity  either  on  the  right 
(S  >  0)  or  on  the  left  (S  <  0),  we  conclude  that  there  is  either  a  root 


25 


i2  >  S  >  0  or  a  root  i  <  -1,  S  <  0.  We  shall  consider  the  case 
for  S  >  0  (the  moment  of  force  tends  to  overturn  the  shell).  The 
other  case  is  quite  similar. 


Since  we  are  interested  in  the  ratio  of  two  integrals,,  we  may- 
neglect  the  factor  of  l/2s.  Then 


Ja 

f 


zn  a  z 


zn  1 


/ 


di 


10.3 


where 


*2  >1  >iQ  >J 


Making  the  substitution 


we  get 


sn  u  =  (i-^1)/(i0-i1) 
k2  = 

™  K 

in  =  (l/K)  J'  (ix  +  sn2u)n  du 


=  (l/K)  i*"J  (i0-ii)J(pA2J 
summing  on  j  and  where 


K 


2j  ,, 

sn  u  du  . 


Now  (see  reference  13) 


10.lt 


10.5 


A0  =  K 

A2  =  k-2  (E  -  K) 

^  AgJ  +  d-g3)  ^ 

(2J+1)  k 


10.6 
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Hence ,  since  there  integrals  are  divided  by  K  to  form  the  averages, 
&  is  a  function  of  iQ,  j8^,  k  ,  and  E/K. 

Finally  we  note  that,  if  we  write  for  the  a’s  in  equation  $.1 
aij  “  Cijn  £ 

summing  on  n  where  the  c's  are  constants, 

"In 


a . .  =  c. .  i 
ij  ijn 


10.7 


and  we  have  the  secular  equations  a  functions  of  H,  0,  and  with 
E/K  as  the  only  transcendental  function. 


One  further  item  is  of  interest  in  connection  with  a  cubic  f .  That 
is  the  geography  of  the  H,  §  plane  (Figure  2).  The  curve  which  bounds 
the  region  of  permissible  motion  has  three  parts:  pure  slow  precession, 
pure  fast  precession  with  p  less  than  90°;  and  pure  fast  precession  with 
p  greater  than  90°.  The  lines  for  one  root  =  1(0°)  and  for  one 
root  -  -1  (lfiO°)  are  drawn  and  are  seen  not  to  intersect  (the  motion 
cannot  oscillate  between  0  and  180°  unless  the  spin  is  zero) .  One  other 
line  is  included.  That  is  the  line  for  H  =  V(i^),  JT  =  Points  on 

this  line  are  those  conditions  which  correspond  to  the  case  of  an 
initially  stationary  shell  axis.  These  are  the  initial  conditions  for 
shell  fired  from  aircraft. 


11.  A  ONE-DIMENSIONAL  EXAMPLE 


The  following  is  given  to  illustrate  the  techniques  described  thus 
far.  Consider  the  one -dimensional  equation 

x”  +  (a  +  b  x2)x*  +  x5  =  0, 

We  need  only  define  one  secular  function,  namely 

H  ■=  (l/2)(x’f  +  (l/4)(x4)  , 

Which  is,  as  usual,  the  sum  of  the  kinetic  and  the  potential  energy. 
Then 

H*  =  -  (a  +  bx2)(H  -  (l/2)x^) 


y 


x1  =  y2H  -  (l/2)x^ 

Define  two  new  variables  r  and  0  by 
x  =  rcos  0 
r^  =  . 
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H,  (D  PLANE 
Figure  2 


Then 


/  2  1  2  2  2 
0s  =  r  sj  (l  +  cos  0)/2  -  (a  +  br  cos  9)(l  +  cos  0)  sin©  cos© 

r*  =  -  (l/2)(a  +  b  r2  cos^Q)(l  -  cos^0)r 


This  is  an  interesting  case.  If  r  is  small  we  cannot  ignore  the 
last  term,  in  the  0*  equation.  Indeed  for  small  r  the  motion  becomes 
oyer-damped.  However ,  we  shall  start  with  a  value  of  r  which  is  large 
enough  to  allow  dropping  the  second  term  in  the  Qf  equation.  This  is 
on  the  assumption  that  a  +  br  is  smaller  than  r /  -J~2 


With  this  as sumption ,  we  let 

cos  0  =  cn(u)  (cn  =  cosine  amplitude) 
in  which  case  the  secularized  r*  equation  becomes 

K 

J2  2  4  g 

(a  +  br  cn  u  x  -  acn  u  ~  ben  u)du 

0 


=  -r  [a/3  +  b(2E/K  -  l)/5  ]  . 

GEDA  solutions  were  run  for  two  specific  cases  (Figures  3  and  4) . 

The  oscillating  curve  is  the  actual  solution.  The  envelopes  are  solutions 
of  the  secular  equation  and  give  the  amplitude  of  the  oscillation.  The 
agreement  is  good. 


12.  THE  MURPHY  NONLINEAR  TREATMENT 
7 

C,  H.  Murphy  has  been  quite  successful  in  applying  the  methods  of 
nonlinear  vibration  theory  to  the  case  of  shell  motion.  It  is  desirable 
to  tie  his  methods  to  the  methods  of  Hamiltonian  mechanics.  To  do  so, 
we  shall  consider  a  somewhat  different  form  of  the  equations  than  that 
of  equations  2.1,  and  for  simplicity  we  shall  consider  nonlinearity  in 
the  Magnus  term  alone. 

We  take  as  the  fundamental  equation  of  motion  (after  Murphy) 

V »  +  (#-  iv)  V  +  ( -M  -  1W)  X  =  0  11.1 

With  v  =  (A/B)v 
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X  +  x=o 

_ H 


I 

J 


I 

j^/X+(X*-|^)X+X*=0 


•'*  “  JL  +  k2  JH 
M  *  k2'2  JM 

T  “  JL  -  V2  JT  *  T1  +  T2 

In  these  equations  we  have  neglected  the  geometric  nonlinearities. 

These  could  be  added,  but  the  algebra  would  become  more  difficult. 

We  now  define  the  variables 

H  -  (1/2)  (X»  I»  -  M  U)  Vv2 

0  =  (U*-  U»)/2iv  -  xI/2  11.2 

These  variables  are  closely  related  to  our  original  definitions.  One 
difference  should  be  noted,  namely,  0"  =0"  , ,-l.  ThiB  means  that  a 

negative  §  is  associated  with  the  slow  precession,  and  a  positive  <j),  with 
fast  precession. 

The  differential  equations  in  terms  of  the  new  variables  are 

+  U/2)  +  TU 


(U)« 

Let 


Then 


and 


H»  =  -2^(H  +  MXX/2v2)  +  t[0  +  XX/2]  11.3 

(sin25)!  =  +  2v  V2sin26  (H  +  sin25/8s)-(0+(sin25)/2)2 
+  Vcr  *\f  (sin^.^  -  sin25)(sin2S-sin28Q)  ' 


(K^  +  Kg)2  =  sin2 8^ 
-  IC,)2  =  sin25Q 


12.4 


sin26  =  K^2  +  Kg2  +  2K^Kg  cos  © 


12.5 


2H  -  (1  +  a)  § 

2 

a 


12.6 
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Using  12.5  and  12,6  in  12.2  and  taking  averages  with  respect  to  0, 


by  using  2 it 


(l/2rt)  J*  (ICj2  +  KpP  +  2¥L^K^  cos  9)  dO  =  K^2  +  K^' 


12.7 


O 


(l/2n)  J'  (Kx2  +  ¥^2  +  21^1^  cos  9)2  d9  =  YL^  +  , 

0 

X,2  +  ?  [T!  +  T2  (*1%  +  2K2a)] 


we  get 

=  2*2 


1,2 


12.8 


4U  -  4,i 


where 


01  2  =  v(2  ±  a)/2 


This  is  Murphy’s  result. 

13.  AIRCRAFT  GUNFIRE 


In  bomber  defense  gunfire  at  high  altitude,  stability  factors  are 
quite  high,  and  the  initial  conditions  are  such  that  there  is  a  large 
slow  precession  and  practically  no  nutation.  Two  things  are  of  interest: 
the  damping  of  the  precession  when  the  nutation  is  zero  and  the  stability 
of  small  nutation.  We  assume  the  stability  factor  is  infinite  which 
allows  us  to  approximate  5.5  and  5-9  by 

$  =  (cos  p)/cos  n  13.1 

H  =  (1/2)  n2 

putting  12.1  Into  3 A  using  the  relation 

cos  6  =  cos  p  cos  n  +  sin  p  sin  n  cos  9  , 
and  keeping  no  terms  in  n  in  the  p’  equation  and  only  linear  terms  in  n 
in  the  n*  equation,  we  get 

p-  =-(JL  -  JT  k^2)  sin  p  13.2 

-2  *«2 

n’  =  -JH  k2  n  -  k^  (n  cos  p  -  sin  p  cos  9) 


The  first  equation  has  been  checked  against  complete  integrations  of 
the  dynamical  problem,  was  found  to  give  good  agreement,  and  is  being  used 
in  the  computation  of  aircraft  firing  tables. 
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To  be  useful  the  second  equation  must  be  secularized  by  averaging 
the  0  term.  The  final  equation  will  be  in  the  form 

2n 

n<  ^  -JHk2 

0 

To  evaluate  the  last  term  let 

JT  -  (cos  p  cos  n  -l-  sin  p  sin  n  cos9)m 

m=0 

z  m  m-  1  „  % 

c^(cos  p  +  m  cos  p  sm  p  cosO  n; 

Then  2it 

(l/2n)  J1  cos  0  dfl  =  ^  m  cosm_1  p  sin  p  15. k 

0 

lk.  NON-RIGID  SHELL 

The  Hamiltonian  formulation  is  well  adapted  to  the  consideration 
of  compound' shell.  In  many  cases  there  are  limitations  imposed  on  some 
of  the  degrees  of  freedom  of  the  component  parts.  We  shall  consider 
one  special  case,  that  of  a  small  mass  which  is  constrained  to  move  in 
a  plane  perpendicular  to  the  shell  axis  and  subject  to  a  linear  re¬ 
storing  force  centered  on  the  shell  axis.  In  linearizing  the  problem 
we  shall  make  extravagant  use  of  canonical  transformations  in  the  hope 
that  the  repetition  will  engender  a  feeling  for  canonical  transformations. 

Consider  the  situation  pictured  in  Figure  5»  The  y  axes  are  space 
fixed  axes.  The  x  axes  are  fixed  in  the  main  shell  body  (the  x^  axis  is 
the  axis  of  symmetry  of  the  main  shell  body) .  The  £  and  axes  define 
the  plane  in  which  the  small  mass  may  move.  The  £  axis  is  taken  in  the 
plane  defined  by  x^,  y  .  The  origin  of  the  £axes  is  at  a  distance  a  from 
the  center  of  gravity  of  the  main  shell  (origin  of  the  x  and  y  axes). 

The  main  shell  body  is  acted  upon  by  a  moment  whose  potential  energy 

is  |i  cos  0.  This  is  again  the  assumption  that  V  is  linear  in  £.  The 

2  2 

spring  has  a  potential  energy  of  (l/2)K((^  +  £2  ) .  The  velocity  com¬ 
ponents  of  the  small  mass  (m)  are  given  by 


-2 


n  -  JT  k2 


n  cos  p  + 


(l/2n) 


sm 


P  ^ 


-2 


f 


cos9  d9 


5k 


14.1 


v1  =  aO  +  -  0  cos  9 

•  • 

Vp  =  f  (a  sin©  +  cos  ©)  +  £2 

«  • 

v5  =  -  0  t;2  sin  ©  -  © 

The  Hamiltonian  for  the  system  is 

H  =  (B/2)(©2  +  f  sln2e)  +  (A/2)  (t  +  0  cos  ©)2  +  \i  cos  © 

14.2 

+  (m/2)  (v^2  +  y22  +  v?2)  +  (k/2)(^2  +  ^2)  , 

where  A  is  the  axial  moment  of  inertia  of  the  main  shell  body  and  B  the 
transverse  moment  of  inertia. 

The  five  momenta  conjugate  to  the  coordinates  £  ,  '£2,  ^  ©,  and  0 
are  respectively: 


•  » 

pt  =  A(t  +  0  cos  ©)  14.3 

«  2  *.  * 

Pe  =  B©  +  ap^  +  m(^  ©  +  m£^£2  sin  ©  0 

d  2  2  2' 

0  =  B  sin  ©0  +  p^cos  ©  +  £2  sin  0  m  0  £-^2m  sin  9 

cos  ©p  +  (a  sin©  +  ^  cos  ©)  p2 

canonical  transformations  with  a 
%ld  ’  ^new^ '  ®le  trails f o rmat ion 


As  a  simplification  in  the  notation  only  those  conjugate  pairs  of  vari¬ 
ables  actually  changed  by  the  transformation  will  be  written  in  the 
generating  fundtion  and  in  the  result.  Those  variables  not  included  will 
be  assumed  to  be  subject  to  the  identity  transformation. 


^2 

We  shall  make  repeated  use  of 
generating  function  of  the  form  F( 
equations  for  this  function  are 


new 


dF 


new 


Pi 


old 


dF 


old 
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For  our  first . transformation  we  take  as  the  generating  function 


It. 5 


P  =  0P3  +  (+  +  0)PU  ,  14. t 

which  gives  us  two  new  coordinates  and  momenta 

q?  =  0 
=  t  +  0 

p3  =  p0  - 

pt  =  p* 

Essentially  this  trans f ormat ion  removes  the  spin  from  the  momentum  p^. 

We  now  define  a  second  transformation  with  generating  function 
F  =  sin  9  sin  p  -  sin  9  cos  q^  Pg  +  (Cq  sin  +  £2  GOS  q3)p7 

+  (-Cjcos  q?  +  ^  sin  %)  p8 

which  gives  the  transformation  equations 
Pq  =  cos  9  sin  q^  p^  -  cos  9  cos  q^  Pg 

p^  =  sin  9  cos  q^  p^  +  sin  9  sin  q^  Pg  +  (f^  cos  q^  -  C2  sin  q^)  Py 

+  sin  q3  +  i2  cos  q?)  Pg 

P1  =  P7  sin  q3  “  p8  cos  q3 

P2  =  cos  q^  +  pQ  sin 


14.6 


it. 7 


q._  =  sin  9  sin  q,. 

5  5 

qg  =  -  sin  9  cos  q^ 
qy  —  sin  q^  +  C g  cos  q^ 
qg  =  cos  q5  +  £2  sin  q? 

This  transformation  introduces  the  direction  cosines  of  the  shell  axis 
(q  and  qg)  and  the  direction  cosines  of  the  vector  between  the  center 
of  force  and  the  small  mass  (q^  and  Qg)*  These  direction  cosines  are 
for  the  y  ,  y2  axes. 


57 


their 


To  simplify  the  situation,  we  shall  consider  q^  g  7  g 

’  2 

derivatives,  and  m  to  be  small.  We  shall  let  K  =  inf  and  assume  that 
the  frequency  of  the  spring  (f)  does  not  vanish  as  m  — *-0 .  Note  that 

this  makes  p„  and  p0  second  order  terms.  Further  we  shall  keep  only 

7  o 

up  to  third  order  terras  in  the  Hamiltonian.  To  do  this  note  that  to 
second  order  terms 


B0  =  (p  -  ap^)  sin  -  (pg  -  apg)  cos  q^ 


14.8 


Bsin  0  f>  =  p. 


sin© 


4  2 

Therefore,  to  third  order  terras 


+  (p  -  ap^)  cos  q^  +  (pg  -  apg)  sin  q^. 


H  =  (1/2B)  J  (p5  ~  ap?)2  +  (p6  -  apg)2  +  p^  ^  -(p5  -  ap^qg  +  (Pg  -  apgjq^] 

>•  O 


+  P 


2  "  2 
2  S  +  % 


2  A  P4 


p  2  2.1  /  2  2v  mf  /  2  2\ 

2  S  +q7  2m  P7  +  P8  2  ^7  +<^8 


For  a  final  transformation  we  shall  use  the  generating  function 
F  =  P9  +  lg  P1Q  +  (aq5  +  q^P-^  +  (aqg  +  qg)pl2 
which  gives  the  transformation  equations 

*9  =  «5 
*io  =  *6 

111  -  aq5  +  q? 

*12  =  aq6  +  qg 
P5  =  P9  +  aPll 
p6  =  p1Q  +  apl£ 

p7  =  P11 


14.9 


14.10 


14.11 


p8 


•12 


This  transformation  changes  the  q7  g  coordinates  from  a  moving  origin 
to  a  fixed  origin  (q  and  g)* 
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Nov  the  Hamiltonian  is 

r  pU  2  tl2 1  P 

H  =  (1/2B)  [(p9  -  q.10  gi)  +  (p1Q  t  %  £)  J  +  (1/2A)  p^ 

-  (n/2)  (^2  +  q102)  +  (l/2m)  (p^  +  p^2) 


(mf2/2)  (< 


L11 


+  (q- 


l12 


aq. 


10 


)2;  . 


If  ve  let  X  =  q^  +  iq^o  and-  z  =  Ppp  +  ^~^±2  an^  no^e  that  =  Av  =  Bv, 
the  Hamilton  equations  give 


X"  -  iv  Is  -  (p/B)X.  =  (amf  /B)  (z  -  al) 


(I  ,  J.2-  j,2, 

z  +  x  z  =  ax  X 


14.13 


Stability  of  this  system  is  equivalent  to  requiring  that  the  roots 
(r)  of  the  determinant  equation 


r2  -  ivr  -  (p/b)  +  (a2rnf2/B) 


af 


amf2/B 

2  j,2 

r  +  f 


=  0 


14.14 


are  pure  imaginary  numbers,  This  vlll  be  the  ease  (for  small  m)  at  least 
If  the  roots  for  m  =  0  are  imaginary  (usual  stability)  and  if  the  roots 
for  m  =  0  are  not  close  together  (resonance).  The  requirement  that  the 
roots  be  pure  imaginary  is  due  to  the  form  of  the  quantlc  polynomial 
represented  by  the  determinant.  If  any  root  is  a  +  ip,  there  Is  also  a 
root  -a  +  ip,  and  one  of  these  roots  introduces  an  unstable  mode. 


While  a  judicious  guess  could  give  equation  14.14,  not  all  cases  of 
restrained  motion  are  as  simple  as  the  preceding.  In  the  more  complicated 
cases,  the  Hamiltonian  method  helps  to  avoid  errors  in  the  formulation  of 
the  equations  of  motion. 


HARRY  L.  REED,  JR, 
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The  exterior  ballistic  motion  of  shell  is  considered  by  the  methods  of 
classical  mechanics.  These  methods  prove  to  be  powerful  tools  for  a  qualitative 
analysis  of  the  non-linear  motion  of  shell  and,  in  certain  special  cases,  for 
obtaining  quantitative  results  with  a  minimum  of  computational  difficulty.  The 
solution  thus  obtained  for  the  aircraft  gunfire  problem  has  proven  to  be  especially 
useful. 
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